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Abstract 

In papers on primary state diffusion (Percival 1994, 1995), numerical estimates sug- 
gested that fluctuations in the space-time metric on the scale of the Planck time 
lO^^'^s) could be detected using atom interferometers. In this paper we flrst 
specify a stochastic metric obtained from fluctuations that propagate with the ve- 
locity of light, and then develop the non-Markovian quantum state diffusion theory 
required to estimate the resulting decoherence effects on a model matter interferom- 
eter. Both commuting and non-commuting fluctuations are considered. The effects 
of the latter are so large that if they applied to some real atom interferometry ex- 
periments they would have suppressed the observed interference. The model is too 
crude to conclude that such fluctuations do not exist, but it does demonstrate that 
the small numerical value of the Planck time does not alone prevent experimental 
access to Planck-scale phenomena in the laboratory. 
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1. Introduction 

Two earlier papers on primary state diffusion, PSDl (Percival 1994), and PSD2 
(1995) describe an alternative quantum theory, and a brief guide to the literature 
in the field. In PSDl the theory had one free parameter, a time tq, which in PSD2 
is close to the Planck time. Thus PSD2 has no free parameters. This paper may 
be considered as a sequel to PSDl and PSD2, in which the proposal to use atom 
interferometry to test the validity of the theory is worked out in greater detail, with 
both commuting and non-commuting fluctuations. There is a large literature on 
experimental constraints on alternative quantum theories, but most are concerned 
with estimates or modifications of free parameters (e.g. Ghirardi, Rimini and Weber 
1986, Pearle and Squires 1994, 1995). Exceptions are the theories of (Diosi 1987, 
1992) and (Penrose 1986, 1996), and of Ellis and his collaborators discussed below. 

We also have a purpose quite separate from any alternative quantum theories, to 
demonstrate that atom interferometers might be used to put an experimental bound 
on Planck scale space-time fluctuations (Ellis et al. 1984). 

It was principally Einstein's 1905 theory of Brownian motion and the subsequent 
experiments of Perrin and others that established without question the reality of 
atoms (Einstein 1956, Pais 1965). Brownian motion is a diffusion process, and 
because of this, measurements on a macroscopic scale could be used to determine 
quantities on an atomic scale. So by analogy with Brownian motion, we should look 
for a diffusion process which enables us to determine quantities on the Planck scale 
by experiments on the atomic scale. Space-time fluctuations produce a diffusion 
in quantum amplitudes (PSD2). We use a relativistic theory of non-commuting 
space-time fluctuations in two dimensions and its detailed application to a simple 
model of an atom interferometer, to show that the small numerical value of the 
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Planck time does not alone prevent experimental access to Planck scale phenomena 
in the laboratory. Section 2 discusses the relation between Brownian motion and 
space-time fluctuations. 

For flnite intervals that are much longer than tq, the space-time fluctuations are 
dominated by drift, and the space-time will appear to be nearly smooth. Neverthe- 
less Planck scale fluctuations contribute to the action integrals in the exponents of 
all quantum path summations. In principle, and possibly in practice, these could 
be detected by atom interferometry. If we could detect the fluctuations we could 
measure the value of tq, which would provide valuable information about dynamics 
on Planck scales. 

In order to do this we treat simple examples of linear non-Markovian quantum 
state difl"usion theory. We treat three cases: commuting delta-fluctuations, commut- 
ing propagating fluctuations and non-commuting propagating fluctuations. Delta- 
fluctuations are introduced by Ghirardi, Grassi and Pearle (1990), whereas non- 
Markovian fluctuations and (different) operator fluctuations are discussed by Pearle 
(1993). 

Quantum state diffusion, or QSD, represents an open quantum system by a pure 
state diffusing in Hilbert space, as a practical alternative to the representation 
in terms of a density operator (Gisin and Percival, 1992 1993a,b, Percival 1994). 
Our formal theory and applications are based on the linear theory of quantum 
state diffusion which was developed independently of the nonlinear theory and for 
a different purpose (Barchielli and Belavkin 1991, Goetsch and Graham 1994). The 
nonlinear theory has the advantages that it is realistic in the sense that a physical 
system is represented directly by a pure quantum state. It can be applied to single 
runs of a laboratory experiment, and is excellent for computations (Schack, Brun 
and Percival, 1995). The linear theory does not have these advantages, but is simpler 
analytically (Strunz 1996), which is why we use it here. Section 3 introduces linear 
Markovian QSD with complex fluctuations in the notation of the usual nonlinear 
theory of Gisin and Percival, and then introduces QSD with non-Markovian simple 
repeated and persistent fluctuations. 

Special relativity was essential to quantum state diffusion as a fundamental theory 
(Gisin 1984), yet it has proved difficult to develop a satisfactory relativistic the- 
ory. Nonrelativistic primary state diffusion (PSDl, PSD2) is one of the realistic 
quantum theories, in which classical systems, quantum systems and the quantum 
measurement process are all represented by state vectors that satisfy the Schrodinger 
equation to a good approximation on atomic scales. They localize on a macroscopic 
scale through a universal diffusion process, with the same mathematical structure as 
nonlinear QSD. In these earlier nonrelativistic formulations this universal diffusion 
is the result of fluctuations which are functions of time only, with a conjecture as 
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to the consequences of relativity. Here we provide simple examples of relativistic 
fluctuations. Any relativistic generalization of a realistic theory like PSD presents 
both conceptual and technical problems (Bell 1990, Diosi 1990, Ghirardi, Grassi & 
Pearle 1990). Wc do not treat the conceptual problems here, but solve some of the 
technical problems that arise because in special relativity the fluctuations depend 
on both space and time. 

According to nonrelativistic primary state diffusion, and by analogy with Brownian 
motion, there is no correlation between the fluctuations at different times. Rela- 
tivistically, to preserve Lorentz covariance, this means that there is no correlation 
over any timelike interval. There is no such direct constraint for spacclikc intervals, 
but nevertheless we shall assume that the only correlations are for null intervals. 
We principally treat the case in which the fluctuations propagate like a scalar fleld 
without interaction, so that there is correlation for all null intervals. For the special 
case of 2D space-time the correlation between any pair of different points on a light 
line is then the same. At the other extreme is the case of delta-correlations in space- 
time (Ghirardi, Grassi and Pearle 1990), where the fluctuations do not propagate 
at all, and which is discussed in the flnal section. 

Among the theories that base a universal decoherence mechanism in quantum me- 
chanics on gravitation is the work of Diosi (1987, 1989 and 1992). He relates funda- 
mental quantum uncertainties in the Newtonian gravitational potential of quantum 
systems (and thus in the time-time component of the metric) to the suppression 
of coherence in their (macroscopic) quantum superpositions. The relevant decoher- 
ence rate is given by the gravitational interaction energy of the superposed states. 
Similar ideas are developed by Penrose (1986, 1996), who relates gravitation and 
wave function collapse. 

The research of Ellis and his collaborators is mainly based on Planck scale physics 
and does not depend on any free parameters. They have suggested that tests of 
quantum gravity can be made using the KK system, and experiments have now 
been made (Ellis, Hagelin, Nanopoulos, & Srednicki 1984, Adler et al, 1995). They 
have also suggested that coherence in SQUIDS should be used (Ellis, Mohanty & 
Nanopoulos 1990). According to PSD the energy difference in both cases is far too 
small for the decoherence due to Planck scale space-time fluctuations to be detected, 
so such experiments might help to distinguish the two theories. Their 1984 paper 
also suggested that neutron interferometry shoiild be used, and the experiments de- 
scribed here are in part a sequel to that suggestion. The 1990 paper emphasized the 
absence of a reliable computational scheme for quantum gravity. This is consistent 
with our approach, which starts from the relatively simple conditions that must 
be satisfied by realistic quantum theories with no free parameters. It seems that 
the problem of uniting quantum gravity with quantum foundations can be attacked 
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from both sides (Amelino-Camelia, Ellis, J. & Mavromatos 1996). 

Section 4 presents a general formal theory with a Dyson expansion for the state 
vectors of the ensemble, which leads to an expansion for the density operator. The 
application of this general theory to particular cases depends strongly on the nature 
of the problem. Section 5 shows how the selfadjoint part of the drift operator is 
determined by the trace condition on the density operator, and the following section 
6 contains expressions for means over products of elementary fluctuations, in which 
the first effects of noncommuting fluctuations are seen. 

Section 7 treats in detail a crude model of a matter interferometer with one space di- 
mension. For this model propagating commuting fluctuations produce no effect, due 
to cancelations between the two arms of the interferometer, whereas non-commuting 
fluctuations and delta-fluctuations lead to a suppression of interference. The spe- 
cial nature of this model leads to some helpful simpliflcations. Section 8 applies 
the theory of the model to the experiment of Kasevich and Chu (1992), and con- 
cludes that it puts severe constraints on the magnitude and nature of the space-time 
fluctuations, subject to the assumptions made. 

2. Brownian motion and space-time fluctuations 

The vertical displacement of a typical Brownian particle is made up of two parts: 
fluctuations and drift. The fluctuations are due to collisions with molecules. These 
are processes on the atomic scale, producing random upward and downward dis- 
placements. The drift, due to the force of the Earth's gravitational field, is uniformly 
downwards. In a time interval At, the mean square displacement M|Asp due to 
the diffusion is proportional to At, so a typical displacement is proportional to the 
square root of the time. In the same time interval the drift is proportional to the 
time itself. Consequently for sufficiently short time intervals the motion is dom- 
inated by diffusion, and for longer intervals it is dominated by the drift. There 
is a characteristic time interval r that marks the approximate boundary between 
the smaller time intervals At for which the diffusion dominates and the larger time 
intervals for which the drift dominates. The measured value of r provided valuable 
information about atoms. For times T much longer than r, the fluctuations are 
typically about {r/Ty/'^ of the displacement due to the drift. 

Space-time is smooth and curved on large scales. The change in the motion of 
particles as a result of the curvature of space-time is seen as the force of gravity. 
The curvature also produces changes in the proper time measured by clocks. But 
space-time cannot be smooth on Planck scales, because of quantum fluctuations. 

For Brownian motion, Einstein was able to obtain relatively simple relations with- 
out detailed dynamics on the atomic scale, which no-one in his time understood 
sufficiently. The corresponding relations for space-time ffuctuations depend on the 
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detailed dynamics of space-time on Planck scales, like the method of quantization, 
and the topology of space-time foam (Hawking 1982, Ellis et al 1984). Ellis and 

his collaborators have based their approach to alternative quantum theories on this 
dynamics. Here, following Einstein, we try to obtain relatively simple relations 
without the detailed dynamics. 

We treat space-time classically, but assume that quantum effects produce stochastic 
fluctuations on the scale of the Planck time. Proper time is then made up of two 
parts: fluctuations and drift. The fluctuations are due to Planck scale quantum 
processes, producing random changes in the proper time. This is added to the drift, 
which is just the proper time when space-time is flat. Consequently for sufliciently 
short time intervals near to the Planck time the proper time fluctuates strongly. 
For longer intervals it is dominated by the drift and looks smooth. There is a time 
To which marks the approximate boundary between the smaller time intervals At 
for which the diffusion dominates and the larger time intervals for which the drift 
dominates. It is generally assumed that the space-time fluctuations are connected 
with quantum gravity, so we would expect tq to be within a few orders of magnitude 
of the Planck time, about 5 x 10~'*'^s. 

For signiflcantly larger scales, a proper time interval As for a timelike segment of 
space-time is represented to a good approximation by 

As = As + ro/'A^(x), (2.1) 

where |A<^(,t)P = As. The barred proper time is for the Minkowski metric, and 
the fluctuations A^(x) depend on the space-time point x. Notice that this is not 
the same as in PSD2, which considers only the time-time components of the metric. 
The formulation here corresponds to the multiplication of a smooth metric by an 
external fluctuating factor. Thus the fluctuations produce a conformal change in an 
otherwise smooth (e.g. Penrose and Rindler, 1984), as proposed by Sanchez-Gomez 
(1994). This leaves null intervals unchanged, and so has the consequence that flelds 
like light in a vacuum, corresponding to particles of zero rest-mass, are unaffected. 

There are many theories of space-time with non-commuting metrics (Connes 1995) 
or an equivalent, including superstring theory (Green, Schwarz and Witten, 1987), 

and supersymmetric GUT (Lopez, Nanopoulos and Zichichi 1994). These are in- 
spired by the need to unite particle theory with gravity. The metrics chosen for this 
paper were not based on any of these, but were obtained by criteria of simplicity 
and accessibility by measurements, from the ideas of alternative quantum theories 
in general and primary state diffusion theory in particular. The non-commuting 
operators are in an isospace, which is attached to the space-time and not to the 
matter in it. As a consequence the factor in the isospace for the density operators 
of matter is the trivial unit operator, which remains unchanged by the dynamics. 
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3. Linear Markovian and non-Markovian QSD 



Linear QSD is obtained from the nonlinear theory by omitting the nonhnear terms. 
In this theory, there is no direct correspondence between the quantum states and the 
states of physical systems, so it is not a realistic theory in the sense used by Einstein 
Podolsky and Rosen (1935) and by Bell (1987). The states are not normalized, and 
each state carries a weight which is given by its norm. However the density operator 
is still given by the same expression as for nonlinear QSD, that is 

p{t)=M\i;{t))m)l (3.1) 

where M represents a mean over the ensemble. This is discussed by Ghirardi, 
Grassi and Pearle (1990). It is important to recognize that the normalized states 
of nonlinear QSD are not given by normalizing the states of linear QSD, but that 
corresponding to every linear QSD theory there is a nonlinear QSD theory which 
gives the same density operator and the same master equation. 

Here we introduce the linear theory in the notation of the earlier Markovian non- 
linear quantum state diffusion theory, without the complications that appear later 
in this paper. First consider the usual case in which each fluctuation is a complex 
scalar d^j{t) which is applied only once, between the times t and t + dt, with the 
ensemble mean orthonormality properties 

Md^j{r)dek{t) = SjkSrtdt 

where by the conventions of the Ito calculus, this is to the lowest order in dt. 

The linear QSD equations are then obtained from the nonlinear equations in Gisin 
and Percival (1993a) by removing the nonlinear parts: 

|dV^) = [- (iA)i/(i)di- i^L](t)L,(t)dt + ^L,(t)dC,(t)]|V^), (3.3) 

3 3 

where the A^j (t) are Ito stochastic differentials. Prom now on all the state diffusion 
equations are linear QSD equations. 

Unlike nonlinear QSD, the linear equations do not preserve the norm of the state 
vector, but nevertheless the density operator obtained from the weighted mean (3.1) 
still satisfies the master equation 



dp{t)/dt = -{i/n)[H{t),p{t)] 

+ E [ - '2L]m,{t)p{t) - ip{t)L]{t)L{t) + Lj{t)p{t)L]{t)_ . 

(3.4) 



7 



The general form of linear quantum state diffusion equation is 

= dG{t)m)) = dR{tmt)) + mt)\m), (3.5) 

where dG{t) is the linear differential evolution operator, made up of a drift part 
dR{t) and a fluctuation part dF{t) whose ensemble mean is zero: 

MdF(t) = 0. (3.6) 



The fluctuation operator dF{t) we assume to be given by the physics of the problem. 
For Markovian QSD it is given by the last sum of equation (3.3). The drift operator 
R{t) is differentiable so dR{t) = R{t)dt, and R{t) can be divided into its self-adjoint 
and skew-adjoint parts, the latter being given by the Hamiltonian: 

R{t) = Rj,{t) + iRi{t) 
Ri{t) = -H{t)/n 

The self- adjoint part of the drift operator is obtained from the so-called trace con- 
dition, that the trace of the density operator p must be conserved for arbitrary p. 
For Markovian QSD this gives 

= Tr dp 

= Tr[|V^)(dV^| + |dV^)(V^| + |dV^)(dV^|] 

(3.8) 



= Tr p\^R + R'^ + J2l]L 

3 

3 



from which the linear Markovian QSD equation (3.3) follows. In this way the trace 
condition is used to obtain both the self-adjoint part of the drift term and the 
resultant master equation. Sections 5 and 7 use the trace condition in the same 
way for non-Markovian noncommuting fluctuations. 

For linear Markovian QSD the differential fluctuation operator dF{t), has the form 

dF{t) = J2Lj{mj{t). (3.9) 

3 

F{t) is not differentiable, despite the Ito notation that suggests that it is. 
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From the stochastic relations (3.2) we obtain for the differential fluctuation opera- 
tors themselves: 

MdF{t) = 0, M{dF{r)dF{t)) = 0, M{dF{r)dF^t)) = X{r)6rtdt, 

(Markov) ''^ 

where X(t) is a time-dependent operator. Equations (3.10) can be taken as the 
definition of Markovian QSD and are particularly useful because they are valid 
when the fluctuations d^ are themselves noncommuting operators, and they can be 
modifled to deflne non- Markovian QSD with repeated or persistent fluctuations. 

Now we remove the Markovian condition, so that the fluctuations can repeat or 
persist over a period of time. We are interested in applications to laboratory exper- 
iments, for which the effects of the fluctuation are small, and so a Dyson expansion 
can be used for the fluctuating quantum states. The measured effects appear in the 
density operator, which is derived by taking the ensemble mean over the perturbed 
states. 

If any fluctuation affects the system at more than one time, then there is correlation 
between the fluctuations at different times, and the diffusion is non-Markovian. In 
the simplest case, which is applicable to the matter interferometer model of this 
paper, each fluctuation appears just twice. This is the example of simple repeated 
fluctuations 

The stochastic relations are then 

MdF{t) = 0, M(dF(i)dF(r)) = 0, 

M{dF{t)dF\r)) = {X'>{r)6tr + X+{r)6t+^t),rX-{r)6t-^t),r)dt, (3.11) 

(simple repeated) 

where a fluctuation which affects the system at time t also affects it at the later time 
t^{t) and the earlier time t~{t), and at no other times. These repeated fluctuations 
are needed for the two-dimensional matter interferometer model of Section 7. 

The theory of simple repeated fluctuations can be generalized to include several 
repetitions and continuously variable time delays. 

4. Dyson expansion 

By standard formal Dyson perturbation methods, if K is the propagator for the 
state vector, so that 

m)) = Km\m), (4.1) 
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then ^ 

K{t,0) =rexp dG{r) 

= rexpG(i) (4.2) 

= ^i^:(-)(i,o), 

n 

where T is the time-ordering operator. 
We assume that for time t = 

m = (4.3) 

for all fluctuations ^, and that 

R{0) = 0. (4.4) 
These are conventions which simplify the analysis. 

In the second line of (4.2) the operator G(t) in the exponent represents the integral 

G{t) = / dG{r) (4.5) 

Jr=0 

and the time ordering operator T refers to the ordering of the operators dG{r). 

When all operators commute, the integrals in the exponent can be reordered in 
any way, and in particular all the operators corresponding to each independent 
fluctuation can be associated with a single time. So for this model the repeated 
QSD equations then become equivalent to ordinary Markovian QSD. A trivial but 
useful example is where G(t) = cl, for all t, in which case the density operator 
remains unchanged. This occurs for the matter interferometry model of Section 7 
with commuting fluctuations. 

On expanding the exponential in (4.2), we get 

K{t,0)=r(l+ f dG{r) + ^ f dG{r) f dG{s) + . . .) 

^ Jr=0 ^- Jr=0 Jr=0 ^ 

= (l+[ dG{r)+ I dG{r) [ dG{s) + . . .) (4.6) 

^ Jr=0 Jr=0 Js=0 ^ 

= (l+ [ dG{r) + / dG(r).G(r) + ...). 

^ Jr=0 Jr=0 ^ 

The initial conditions and iteration for K^'^\t, 0) are given by 

K^'^\t,0)=I, K^^\t,0)=G{t), K^''\t,0)= f dG{r)K^''-^\r,0), 

Jo 

(4.8) 
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where the limits of integrals are all labeled by the time, so for = 0), 
we have 

^(0) _ ^(1) _ ^(2) _ r dG(r).G(r). (4.9) 

^0 

Despite appearances, K^'^^ cannot usually be integrated to give G{t)^/2, because 
dG{r) does not usually commute with G{r). 

The first few terms in the expansion of the density operator pt at time t are 

Pt = MK{t,0)poK{0,t) 
= M(po 

+ K('^Po+PoK^'^^ + K^'^PoK^'^^ 
+ K^^^po + poK^^^'' + K^^^PoK^^^'' 

+ ...), (K^^) =K(^)(t,0)). 

Section 6 shows that the fluctuations contribute to the diagonal terms K^'^'' pqK^'^^^ 
only. 

5. Trace condition 

The operator G has to be expressed in terms of the drift and fluctuation operators 
R and F, but the the drift terms are not known until the trace condition is applied. 
The only terms in the expansion that can be obtained directly are the fluctuation 
and Hamiltonian terms. The self-adjoint part of the drift operator has to be eval- 
uated simultaneously with the term by term evaluation of the density operator, 
which makes things complicated. 

However it often happens that the only significant term is the first non-zero term. 
This is certainly true for the applications that we have in mind, for which any 
nonzero effect of the fiuctuations would be important. In that case the procedure 
is relatively simple. 

The propagator with the drift terms set to zero is the pure fiuctuation propagator, 
denoted KF{t, 0), for which the entire above theory applies with G replaced by F. 
Let pFt denote the density operator obtained from the fiuctuation propagator alone, 

PFt = MKpit, 0)poKUt, 0). (5.1) 

The expansion of this partial density operator can be simplified, because the off- 
diagonal terms MK^^ pqK^ are zero for different n and n'. This follows from 
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(4.10) 



the theory of unbalanced means presented in the next Section. The expansion then 
becomes 

PFt = Y.MK^^^PoK^^^\ (5.2) 

n 

This is expanded up through the first non-zero term, labeled by n = ni, which is 
then used to obtain the corresponding drift to the same order, using the condition 
that the trace of p is 1 for arbitrary initial po, giving 

= T¥(pi - po) ~ T¥(5po + PoS"^ + MK^p^poKp^^), (5.3) 
so that the self-adjoint part of S is 

^R^-ii^?^^^4"^^ (5.4) 

It follows that any term in the expansion of ppt which is proportional to po has 
no effect on pt because it is canceled by the corresponding drift term which comes 
from the trace condition. This helps to simplify the analysis. 

6. Commuting and non-commuting fluctuations 

In this section wc obtain means over prodiicts of commuting and noncommiiting 
fluctuations. Here we use differences instead of differentials, and then take the small 
time limit, because the limiting processes are subtle. Before taking the limit the 
equalities are correct only to leading order in powers of At. The discrete times t 
are separated by multiples of the interval At. We later use the limit 'lim' which is 
always to be understood in the sense of 'limAt^o'- 

In earlier versions of QSD each fluctuation was supposed to be a complex scalar 
with distribution invariant under a complex rotation around the origin, represented 
by a multiplying factor of modulus unity. A point in the space of such fluctua- 
tions can be represented by one complex parameter. This is generalized to non- 
commuting fluctuations which have the corresponding statistical properties with 
the complex conjugate replaced by an adjoint. The number of complex parameters 
is N and they can be chosen so that the distribution of the fluctuations is invari- 
ant under rotation about the origin in a complex Euclidean parameter space, an 
iso-space independent of space-time. It is helpful to think of the distribution as a 
Gaussian in this parameter space. The means over products, which appear in the 
density operator, always reduce to a real number times the unit operator in the 
iso-space. The operators X of Section 3 have the same trivial factor. 
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If j labels a set of independent fluctuations, then they satisfy the basic stochastic 
equalities 

MA^j{t) = (a) 
MACj{r)Aa{t) = (b) (6.1) 

MA^j{r)A^l{t) = SjkSrtAt (c), 

which should be compared with equation (3.2) for commuting fluctuations. Notice 
that (a) and (b) follow from rotational invariance, since almost every combined 
rotation in the parameter spaces of the fluctuations changes the value of the mean 
unless that value is 0. The same goes for the zero off-diagonal elements in (c), for 
which the diagonal element provides a normalization condition. 

In the same way the rotational invariance shows that an unbalanced mean is zero: 

MACj{t)''AC}{t)"' = (n,m different). (6.2) 

This result is used in eqn. (5.2) to simplify the expansion of Kpit) by removing the 
off-diagonal terms with different n, n' from the expansion (4.10). 

The second order expansion of Kp leads to products of four fluctuations for pp. 
This requires the mean 



MA^(r')A^(r)A^t(i)A^t(i')^ (6.3) 



which is zero unless 



r = t and r' = t' (direct) 
or r = t' and r' = t (exchange). 

These are called the direct and exchange terms as shown. The direct term is in- 
dependent of the commutation properties of the fluctuations. It can be evaluated 
directly to give At^ by first taking the mean over the fluctuations A^{t) and then 
over A^{t'). The exchange term depends on the commutation properties of the 
fluctuations. If they commute, then it is the same as the direct term. Otherwise we 
deal with each case separately. 

In this paper we restrict the detailed theory to the special case oi N = 3 where 
the fluctuations are derived from Pauli matrices cr^. These will be called Pauli 
fluctuations. Then in order to satisfy the basic relations (6.1) the A^{t) are given 
by 



where A^,(r)A^J(t) = SijSrtAt 



(6.5) 
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so that the 3 complex components A$^i{t) are statistically independent fluctuations. 
They can be considered as components of a complex 3-vector fluctuation in an 
iso-space. 

The flrst product in the Dyson expansion that depends on the commutation prop- 
erties is the second order exchange term, which is 

^ (6.6) 

= i(-3)(At)2 = -l(At)2, 

where sums are over all suffixes and arguments. In the last sum there are 3 products 
in which i' = i and af = 1, and 6 products in which i and i' are not equal, with the 
value (cri(j2)^ = — o-| = —1. 

For arbitrary fluctuations that satisfy the conditions of this section, 

MA^(tOA^(t)A^t(^')A^t(^) ^ rj{Atf, (6.7) 

for some constant ij which depends on the commutation properties of the fluctua- 
tions. For commuting fluctuations and for the Pauli fluctuations above we have 

77(commuting) = 1, 77(Pauli) = —1/3. (6-8) 

We do not consider cases with all times identical, for example 

MAmm)^^Ht)^^Ht), (6.9) 

because althoiigh they are of the same order as the direct and exchange terms, their 
contribution tends to zero with At. 

Summarizing for the second order: 

MAar')AC{r)AC\t)AC\t') = {5r't'Srt + r)5rt'5r't){Atf . (6.10) 
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7. Interferometer model 



Here we consider a specific model which is a crude representation of a matter inter- 
ferometer in a fluctuating field. The wave packet in each arm of the interferometer 
is represented by a single state, with label 1 on the left and 2 on the right. These 
states have orthogonal projection operators Pi and P2. For the purposes of this 
paper the size of the wave packets is supposed sufficiently small that they can be 
considered to be at points in space. The theory is worked out for 1 space dimension 
only. The remaining dimensions complicate the theory, but they do not affect the 
orders of magnitude of the results, which is all that we are concerned with here. 
For the atom interferometers of interest, the velocities are typically of order lms~^, 
so we neglect the effects of time dilation due to these velocities here. 

Let the paths of the wave packets be Xj{t),j — 1,2. Choose the origin of time when 
the wave packets separate, and let T be the drift time before they recombine to 
produce an interference pattern, so that 

xi(0)=X2(0), xi{T)=X2{T) = 0. (7.1) 

An interaction representation is used in which the basis for each wave packet is the 
unperturbed wave packet, so that the interaction Hamiltonian is zero. 

We use time units to measure distances, with the velocity of light c = 1, so that 
lns~ 0.3m. For an atom interferometer whose wave packet separations are produced 
by photon recoil, wave-packet separations might be about 0.1ns, whereas the drift 
time T is of the order of Is. 

The delta-fluctuations for the two wave packets are independent of one another, so 
the QSD equation and the master equations are (3.3) and (3.4) with H = 0, Li = 
r^/^Pi, 1/2 = r^/^P2- As a result, the off-diagonal elements of p decay exponentially 
with decay constant F, where F is the inverse of the decoherence time defined in 
PSDl, 

F = Tp-i = {mc^fToh-^ (7.2) 

Thus, for delta-fiuctuations, the interference pattern is suppressed significantly if 
the drift time T of the wave packets is not less than F~^. 

In the case of the commuting propagating fluctuations, for each time t there are 
two impulsive fiuctuations A^+(t) and (t) that reach wave packet 1 at time 
t, which last for a time At and which are labeled by their time of arrival at wave 
packet 1. They propagate with the velocity of light to the right from 1 to 2 for 
A^+{t) and to the left from 2 to 1 for A(-{t) . 

In the following equations we approximate using the small velocities of the wave- 
packets relative to light. It follows that to a good approximation the distance 
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between the wave packets can be assumed constant for the time of propagation of 
hght between them, or for a small multiple of that time. 

If x{t) is the distance between the wave packets at time t (in time units) and 
tj is the time that a + fluctuation passes j, then in this approximation the time 
of propagation can be taken as x{t)^ where t is ti or t2 or any time in between. 
Similarly in the propagation of a + fluctuation from 1 to 2, followed immediately 
by the reverse propagation of a - fluctuation from 2 to 1, takes time 2x(t), where t 
is any time during the propagation of either fluctuation between the wave-packets. 
In the following we use the deflnitions 

t+=t + x{t), t~=t-x{t), t++=t + 2x{t), t~~ =t-2x{t). 

(7.3) 

It follows that in this approximation 

t+<t'^t<t'~, t+ <t'~ ^t<t'~~, (etc.) (7.4) 

A characteristic function x(condition) is 1 when the condition is satisfled and 
when it is not. The above approximation will be used in the characteristic func- 
tions of eqns (7.9) and (7.10) for the second order expansion with noncommuting 
fluctuations. 

For each time t there are two fluctuations A^'''(t) and A^~{t) that propagate from 
1 to 2 and from 2 to 1 respectively with the velocity of light. The Hamiltonian is 
zero, each fluctuation is applied jiist twice, once on each wave packet. The evolution 
of this model system is represented by a differential fluctuation operator 

AF{t) = [Pi(A^+(t) + + P2(A^+(t-) + Ar (t+)], (7.5) 



The fluctuations A^^ may or may not commute, but even when they do not com- 
mute with each other, they commute with the projection operators and the density 
operator, so they live in an 'iso-space' that is distinct from ordinary position space. 

When the fluctuations commute it is convenient to deflne a total fluctuation operator 
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at time T' given by 



T' 







T' 



= lim^[Pi(A^+(t) + ACm + P2(A^+(i) + AC{t)] 





T' 



(7.6) 



= /lim^(Ae+(i) + Ar(i)) 








where / is the unit operator. 

For commuting fluctuations we can take the exponential without time ordering to 
get the propagator Kp, which is 



and proportional to the unit operator, so by Section 5 the full propagator is unaf- 
fected by the fluctuations. So the density operator remains at po in the interaction 
representation, and the interference pattern is unchanged when T' = T. 

Now consider the perturbation expansion (4.10) of Pt for non-commuting fluctua- 
tions. The commutation properties do not affect the first two terms in the expansion, 
which therefore make zero change in the density operator. 

(2) 

The first significant term in the expansion is therefore Pprp, which is obtained using 
(6.10), and the definitions 




(7.7) 



Pdg = PipoPi + P2P0P2 



pod = P1P0P2 + P2poPl- 



(7.8) 
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It is 



= lim Mx{r < r')SF{r')SF{r)pox{t < t')SF'' {t)SF'' {t') 

rr'W 

= lim {At)^x{r < r')x{t < t') 4pdg(5rt^r't' + vSrt'Sr't) 

rr'W 

+ Pod[{Sr,t- + Sr-,t){Sr',t'- +^r'-,t') 

+ Vi^r',t- + ^r'-,t)i^r,t'- + ^r",*')] 

= \imJ2{^tfx{t<t') [4pdg[l + 0] 



+ Pod[x{t+<t'^)+x{t+<t' ) + x{t-<t'^)) + x{t-<t' ) 

+ rj(x{t'^<t+) + x{t'^<t-) + X{t'-<t+) + X{t'-<t-))] 



= \unJ2{^tfx{t<t') 



w 



4pdg 



+ Pod [xit<t') + xit<t' ) + x{t<t'^^)) + xit<t') 
+ v(x{t'<t) + x{t'^^<t) + x{t'~~<t) + x{t'<t))] 



(7.9) 



The product of characteristic functions with incompatible arguments is zero, and the 
product of a characteristic function with a stronger condition and a characteristic 
function with a weaker condition is equal to the stronger, so 

w 

+ Pod[Mt<t') + X{t<t'~~) + r]x{t'~~<t<t')] 

= \imJ2{Atf[^PoX{t<t')+Pod{ri - l)x{t'~~<t<t')] (7.10) 
w 

dt' I dMpo+ / di' / dipod(r7-l) 

^0 Jo Jt'-2x{t') 

= 2TVo + '2A{r] - l){PiPoP2 + P2P0P1). 

where A is the area enclosed by the two paths of the interferometer in space-time, 
measured in units of time^. The first term is proportional to the initial density 
operator, and so is canceled by the corresponding term in Sr, as shown in Section 
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5. The second term does not affect Trpy, and so does not contribute to Sr. So we 
have 

PT ~ PO + rM(77 - l)(PlP0^2 + P2P0P1) 

PT ~ Po - r^A^{PipoP2 + P2P0P1) (Pauh). ^^'^^^ 

The state diffusion due to the fluctuations tends to suppress the off-diagonal ele- 
ments of the density operator, and thus decoheres the wave-packets in the two arms 
of the interferometer. This will be detectable by a suppression of the interference 
pattern if r^^(l — ry) is sufficiently large. Any value comparable to 1 would be 
enough. 

8. Numerical values and discussion 

At present there is no known positive evidence from experiment that space-time 
fluctuations of any kind suppress the interference of matter interferometers. The 
negative evidence, that any suppression is below the experimental limits, puts a very 
provisional upper bound on the possible value of the fundamental time constant 
To for the space-time fluctuations. The bound is provisional because it depends 
on many assumptions. It is assumed that the two-dimensional model adequately 
represents a real interferometer. 

Of the three examples in the introduction, the commuting fluctuations that prop- 
agate with the velocity of light produce no decoherence when time dilation is ne- 
glected. For the other two, suppose that the experiments put an upper bound of 
10% on the reduction of the off-diagonal elements of the density operator p by the 
space-time fluctuations of eqn (7.11). 

Under these assumptions 



FT < 0.1 (delta-fluctuations) 
r^8A/S < 0.1 (noncommuting fluctuations) 



and so from the value of the decoherence rate F given by (7.2), the time constant 
of the space-time fluctuations is bounded by 



1/ h \2 lo^'^^s'' 

^° ^ ^ ) ~ 2A^T (delta-fluctuations) 

/0.3\i/2/ h \2 10-49g2 
^° ^ ySAJ \A^ ) ~ IM172 (noncommutmg fluctuations), 

(8.2) 

where A is the atomic number and u is the atomic mass unit. 

The best upper bound is proportional to the inverse square of the atomic number A, 
as given in PSDl and PSD2. But the dependence on area and time is different from 
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the formulae in those papers, which is why our conclusions are different. Another 
new feature is that according to the theory presented here, in which the time dilation 
for motion of the atoms is neglected, the propagating fluctuations must not commute 
if there is to be any efi^'ective bound. 

The best bounds are the smallest. They are given by the more massive particles, 
which suggests atom interferometry rather than neutron interferometry, and a rel- 
atively large value of the area A. These conditions were met by the experiment of 
Kasevich and Chu (1992), for which the atom was sodium with A = 23 and the 
area was approximately 10~^^s^ and the time T ^ 50ms so that 

To < 1.8 X lO^^^s (delta-fluctuations) 

tq < 1.8 X 10 s (noncommuting fluctuations). 



These are significantly less than the Planck time of about 5 x 10^ s, which shows 
that this experiment puts a severe bound on possible space-time fluctuations, given 
the assumptions. 

What about the assumptions? Theories of space-time with non-commuting metrics 
or an equivalent include superstring theory, GUT and the theory of Connes. They 
are inspired by the need to unite quantum theory with gravity. The metrics chosen 
for this paper were obtained by the conditions of simplicity and accessibility by 
measurement, and from a perceived need for a realistic quantum theory in the sense 
of Einstein, Podolsky and Rosen (1935) and of Bell (1987). Any similarity to other 
non-commuting metrics is a bonus. 

The main weakness of this paper is the crudity of the model matter interferometer, 
because we have taken the Hamiltonian to be the mass together with projectors onto 
the wave packets in the arms of the interferometer, neglecting the kinetic energy 
terms. This is a non-trivial assumption. 

However, despite these assumptions, we have shown by means of the model that 
Einstein's method of accessing small quantities using diffusion processes, which he 
applied so successfully to the atomic scale using Brownian motion, might also be 
applied today to carry out experiments on Planck scale space-time fluctuations. 
The small value of the Planck time is not an impassable barrier to the precision 
of modern matter interferometers. Even if the experiments were to show that such 
fliictuations are not significant on a Planck scale and below, which looks possible, 
this would be a contribution to a field that has lacked such experimental evidence 
in the past. 
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